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Abstract 
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t^. ' The spectrum of orbitally excited D s mesons is computed in the continuum limit of quenched 

o' ... 

CH ■ lattice QCD. The results are consistent with the interpretation that the narrow resonance in the 
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D s tt° channel discovered by the BABAR Collaboration is a J p = 0+ cs meson. Furthermore, 
within statistical errors, the 1 + — 1~ and the + — 0~ mass splittings are equal, in agreement with 
the chiral multiplet structure predicted by heavy hadron chiral effective theory. On our coarsest 
lattice we present results from the first study of orbitally excited D s mesons with two flavors 
of dynamical quarks, with mass slightly larger than the strange quark mass. These results are 
consistent with the quenched data. 
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INTRODUCTION 



The recent discovery by the BABAR collaboration of a new resonance, with a mass 
around 2.32 GeV and a narrow width, in the Dfir final state [l| has provoked a great deal 
of interest from experimenters and theorists alike. The CLEO collaboration |2| has confirmed 
this resonance. Both experiments interpret this as the lowest lying of the four P-wave states 
the 3 Pq with J p = + . A variety of theory papers have been published 
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12| either supporting this interpretation or presenting alternative hypotheses of "exotic" 
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states. These alternatives are motivated by potential model results |13l. Il4j[. which suggest 
that the cs scalar meson mass is around 2.48 GeV and hence above the DK threshold. One 
could argue the merits of a particular model, but this becomes irrelevant if the spectrum 
can directly determined from QCD. A mass below the DK threshold would also explain the 
narrowness of the state. 

There have been a number of previous lattice QCD calculations of the L = 1 states in 

n 

the D s meson spectrum. Hein et al. [15j obtained 500(80) MeV for the D* — D s mass 
splitting (from their figure 27) using NRQCD at a fixed lattice spacing of 0.18 fm. Lewis 
and Woloshyn |r| obtained 530(15)(5) MeV at a fixed lattice spacing of 0.11 fm also using 
NRQCD. With quenched relativistic heavy quarks, Boyle 17[ obtained a mass splitting 
of 544(20) MeV at a fixed lattice spacing of 0.07 fm. Although the calculation used the 
same ensemble of gauge configurations as this work, the propagators were computed using a 
slightly different action and different definition of the lattice quark mass. Recently Bali |10| 
presented results in the static limit for the heavy quark and obtained a value of 468(43) (24) 
MeV for the scalar-pseudoscalar mass splitting in this limit. 

All the previous lattice QCD calculations of the D s spectrum were done at fixed lattice 
spacing in quenched QCD. In this calculation we take the continuum limit in quenched QCD, 
so that lattice artifacts are under control. The lattice volumes are large enough (greater 
than 1.5 3 fm 3 ) that finite size effects should be small. We also report results from the first 
unquenched calculation at fixed lattice spacing. The remaining systematic uncertainty is 
due to dynamical u and d quarks having unphysically large masses. 

In the heavy quark limit, the spin of the heavy quark decouples from the rest of the 
system. The observable states can be labelled by the total angular momentum of the light 
quark j [18]. The P-wave states have L = 1 and thus j = {|, §}. This combined with the 
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spin of the heavy quark produces two doublets. The j ' — | doublet contains a J = 2 and a 
J = 1 state, the j = \ doublet contains a J = 1 and J = state. The two J = 1 states do 
not have definite charge conjugation and so can mix. On the lattice only the lightest state 
in this channel can be determined easily. 

In the double limit of heavy quark and chiral symmetry, the two heavy light multiplets, 
{0~, l - } and {0 + , 1 + }, are degenerate. The effect of spontaneous chiral symmetry breaking 
is to split these parity partners, such that the mass splittings 1 + — 1~ and + — 0~ are 
equal [2|. This is confirmed by the CLEO Collaboration 2] who obtain splittings of 351(2) 
and 350(1) MeV respectively. It is interesting to explore the extent to which QCD reproduces 
this remarkable agreement. 



LATTICE DETAILS 



The spectrum of cs mesons has been determined on four ensembles of gauge configura- 
tions. Three have different lattice spacings (a) and were generated in the quenched approxi- 
mation, which enables the continuum limit to be taken. The fourth ensemble was generated 
with two degenerate flavours of dynamical quarks, and the lattice spacing matched to that 
of the coarsest quenched ensemble. For all the ensembles, the Wilson gauge action and the 
non-perturbatively 0(a) improved Wilson fermion action were used. The lattice parameters 
are detailed in Table H] The procedures for generatin g th e dynamical ensemble and matching 
to the coarsest quenched ensemble are described in 19j . 

Meson correlation functions were computed with several different heavy quark masses 
which span the charm quark mass, and several light quark masses around the strange quark 
mass. For the dynamical ensemble only one sea quark mass is used, for which the ratio 
mps/mv = 0.70(1) when m sea = "^valence- This corresponds to QCD with two dynamical 
flavours of mass slightly above the strange quark mass. The details of extracting the spec- 
trum from lattice correlation functions, and the results for the S'-wave D meson spectrum 

20|. To measure the mass splittings, the 



for the finer two lattice spacings can be found in 
ratio of correlation functions at large times is fitted. We have checked that the results from 
computing the 1 + — + mass splitting directly is the same, with the same statistical errors as 
obtained by combining the results obtained for the 1 + — 1™, + — - and 1~ — 0~ splittings. 
In quenched QCD, or in simulations with unphysical heavy sea quarks, there is an am- 
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biguity in the determination of the lattice spacing in physical units. For typical simulation 
parameters used today this is estimated to be of the order of 10%. The scale is set through- 



out this calculation from the static quark potential using tq |2l|, |22j. The value of r /a is 
unambiguous for each ensemble, and so is a good choice for comparing results from differ- 
ent ensembles. However, there is no agreed experimental value for r . Sommer originally 
advocated rp = 0.5 fm. The lattice spacing obtained from the K*/K mass ratio (method 
of planes) |23( corresponds to r$ ~ 0.55 fm on the ensembles used in this work. Determi- 
nations of the lattice spacing from the kaon decay constant, the nucleon mass or the rho 



mass correspond to r values ranging from approximately 0.5 to 0.55 fm [19, 20, 24, 25, 26]. 
For these reasons we take the value of r to be 0.55 fm. The analysis has been repeated 
using r = 0.5 fm throughout, and the difference is taken as an estimate of the systematic 
uncertainty in the scale. 

The strange quark mass is set by from the light-light pseudoscalar mass with the experi- 
mental kaon mass as input 21\. Similarly the charm quark mass is set from the heavy-light 
pseudoscalar mass with the experimental D s meson mass as input. 

RESULTS 

For each of the four ensembles, the mass splittings 1 + — 1~ and + — 0~ are equal within 
statistical errors. These results and the continuum extrapolation which is linear in a 2 are 
shown in figure ^ Furthermore, these splittings are also equal in the continuum limit, in 
agreement with heavy hadron chiral effective theory and experiment. At the coarsest lattice 
spacing the effect of introducing sea quarks with a mass close to the strange is to slightly 
lower the 1 + — 1~ and the + — - mass splittings, but this is not statistically significant. 

Shown in figure El and Table |H] is the comparison of the lattice results along with the 
experimentally measured spectrum. The dynamical results appear to have smaller error 
bars and to be systematically higher than the quenched result. These effects are due to 
extrapolating the quenched results and are absent at fixed lattice spacing. The computed 
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1~ — 0~ splitting is too small, a well known failing of the quenched approximation |2£ 
Evidently the dynamical sea quark mass is too large to chan ge t his, as has been observed 
before in the light hadron spectrum on the same ensemble [19j. The lattice results for 
the cs + and the lightest 1 + mesons are consistent, albeit within large statistical and 
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systematic uncertainties, with the masses of the states discovered recently by BABAR and 
CLEO. Although u and d sea-quark effects are not yet properly included, these lattice results 
provide the most reliable computation of the cs spectrum to date. Our errors are too large 
to exclude exotic states based on potential models. However, there is no evidence from 
our lattice QCD calculations that exotics are required to explain the BABAR and CLEO 
discoveries. 

The lattice data was generated on the Cray T3D and T3E systems at 
EPCC supported by, EPSRC grant GR/K41663, PPARC grants GR/L22744 and 
PPA/G/S/1998/00777. We are grateful to the ULgrid project of the University 
of Liverpool for computer time. The authors acknowledge support from EU grant 
HPRN-CT-2000-00145 Hadrons/LatticeQCD, and PPARC grants PPA/G/O/2000/00456, 
PPA/P/S/1998/00255(CMM), PPA/N/S/2000/00217(RDK). 
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FIG. 1: Continuum extrapolation of quenched results (open symbols) and Nf = 2 results in 
dimensionless units at fixed lattice spacing (solid symbols, offset for clarity ) for mass splittings in 
the D s system. 



TABLE II: Comparison of lattice results with experiment. Mass splittings from the -Dif (1969) 
Mass in MeV. CL denotes continuum limit, tq = 0.55 fm is used unless otherwise noted. 



J p 




experiment 


N f = CL 
ro = 0.5 fm 


Nf = CL 


N f = 2 
a' 1 ~ 1.7 GeV 


N f = 
a' 1 ~ 1.7 GeV 
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FIG. 2: Comparison of experimental results with the lattice determinations from the quenched 
continuum limit (open symbols) and the Nf = 2 data at fixed lattice spacing (solid symbols) 
with ro = 0.55 fm. For comparison we also show the quenched result at the same lattice spacing 
(open diamonds). Zero on the vertical scale is set by the 7J^(1969) mass. Also plotted are the 
experimental DK and D*K thresholds. 
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